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Abstract 

This paper presents a new method for analyzing near-field Faraday probe data obtained from Hall 
thrusters. Traditional methods spawned from far-field Faraday probe analysis rely on assumptions that are 
not applicable to near-field Faraday probe data. In particular, arbitrary choices for the point of origin and 
limits of integration have made interpretation of the results difficult. The new method, called iterative 
pathfinding, uses the evolution of the near-field plume with distance to provide feedback for determining 
the location of the point of origin. Although still susceptible to the choice of integration limits, this 
method presents a systematic approach to determining the origin point for calculating the divergence 
angle. The iterative pathfinding method is applied to near-field Faraday probe data taken in a previous 
study from the NASA-300M and NASA-457Mv2 Hall thrusters. Since these two thrusters use centrally 
mounted cathodes, the current density associated with the cathode plume is removed before applying 
iterative pathfinding. A procedure is presented for removing the cathode plume. The results of the 
analysis are compared to far-field probe analysis results. This paper ends with checks on the validity of 
the new method and discussions on the implications of the results. 

Nomenclature 


GRC 

Glenn Research Center 

300M 

NASA-300M 

457Mv2 

NASA-457Mv2 

HARP 

High-speed Axial Reciprocating Probe 

CEX 

Charge Exchange 

SEE 

Secondary Electron Emission 

CD 

Current Density 


Introduction 

Recently, Hall thrusters have been seeing a rise in popularity as primary spacecraft propulsion in the 
United States. Examples include the use of Busek’s 200-W Hall thruster on TacSat-2, (Ref. 1) the use of 
Aerojet’s BPT-4000 Hall thrusters (Ref. 2) on the Air Force Advanced Extremely High Frequency 
satellite program, and NASA mission concepts that use high power Hall thrusters to reduce the cost of a 
manned trip to Mars (Refs. 3 and 4). While the community understands the physics and engineering of the 
Hall thruster well enough to develop flight units, there are still many gaps in that understanding. One of 
those knowledge gaps relates to the transition of Hall thruster plasma plume from the well-defined 
annular jet at the thruster exit to the diffused plume in the far-field. 

To improve understanding of plume transition in a Hall thruster, one approach is to insert electrostatic 
probes into the plume and interrogate for local plasma properties. However, the plasma in the near-field 
of a Hall thruster is highly energetic and will cause stationary probes to overheat, possibly to the point of 
probe failure. In 2000, Haas, et al., pioneered the use of a high-speed axial reciprocating probe (HARP) 
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system to overcome the heating issue by rapidly inserting and retracting probes into and out of the near- 
field region of a Hall thruster (Ref 5). Since then, many new tests involving various electrostatic probes 
have been performed using the same concepts and methods (Refs. 6 to 12). This paper will focus 
specifically on the development of the near-field Faraday probe. A Faraday probe is a flat metal disk, 
sometimes with a guard ring, that is biased well into ion saturation regime of the local plasma in order to 
measure the ion current density. Analysis of the ion current density as a function of position can yield 
information on total ion beam current and plume divergence. 

Currently, the biggest obstacle associated with using the near- field Faraday probe is the ad-hoc nature 
of the existing data analysis methods. There are two outstanding issues with existing analysis methods. 
One, existing methods typically start with assumptions about the point of origin from which the near- field 
plasma emanates. These assumptions are not grounded in physics or plasma properties but on thruster 
geometry. Two, the analysis results are often strongly tied to the integration limits, which are necessary 
due to the otherwise infinite nature of the integration domain. 

A new data analysis method, called iterative pathfinding, is presented in this paper. This new method 
is an attempt to address the first of the aforementioned issues. The general principle of the iterative 
pathfinding method can be summarized in two steps. Step one, using a guess on the point of origin and 
the divergence equation, the direction and spread in the plasma plume is calculated. Step two, data on the 
evolution of the direction and spread in the plume with distance is fed into an algorithm to find a new 
guess on the point of origin. Steps one and two are repeated until convergence is achieved. In order to 
establish the error metrics for evaluating the next guess on the location of the point of origin, an 
assumption on how the near-field plasma evolves with distance is needed. For this study, the near-field 
plasma plume is assumed to undergo free-expansion, with no self-interaction, after exiting the channels of 
the Hall thruster. No assumption is made on the angular distribution of the plasma plume. While the new 
method is still sensitive to the choice of integration limits, the resulting origin points are physics- and 
plasma property-based. Validity of the new method is evaluated by assessing the sensitivity of the results 
to integration limits and comparing said results to those obtained from a far-fleld Faraday probe. 

Past Work 

In 2000, Haas, et al., pioneered the use of electrostatic probes in the near- field of the Hall thruster by 
developing a technical solution to the problem of near- field probe heating (Ref. 5). Shortly after, in 2003, 
Hofer and Gallimore used this technical solution with a near- field Faraday probe to map the near-field ion 
current density of the NASA-173Mv2 (Ref. 7). Being the first near-field Faraday probe study, the authors 
did not have time to fully explore multiple operating conditions and bias voltages. No divergence analysis 
was performed and total beam current was calculated at one axial distance. 

After the pioneering days, near-field Faraday probe studies did not surface again in the literature until 
2008, when Jameson (Ref. 9) and Reid (Ref. 8) independently applied it to two separate H6 Hall thrusters. 
Jameson used the near-field Faraday probe data to calculate total beam current at one axial distance and 
relied on far-field Faraday probe data for divergence calculation. Reid proposed a method for measuring 
the spread in the plume so that comparison can be made across different operating conditions. Reid also 
tried several different schemes for determining the integration limits. 

Since then, three more near- field Faraday probe studies have been published including the 2011 work 
by Shastry on an H6 Hall thruster (Ref. 13), the 2012 work by Hofer, et al., on the H6MS (Ref. 14), and 
the 2012 work by Huang, et al., on the NASA-300M and NASA-457Mv2 (Ref. 1 1). The 2011 work by 
Shastry followed the methodology outlined in Reid’s 2008 work. The 2012 work by Hofer, et al., was 
focused on the application of magnetic shielding to the H6 Hall thruster and did not describe the near- 
field Faraday probe analysis in detail. 

The 2012 work by Huang, et al., on the NASA-300M and NASA-457Mv2 proposed a new physics- 
based data analysis method that borrowed a few concepts from optical ray tracing (Ref. 11). The authors 
noted that the near-field current density profiles in the vicinity of the discharge channel have Gaussian- 
like distributions. They assumed that the near-field plasma plume can be approximated as free-expanding 


NAS A/TM— 20 13-216567 


2 



Gaussian jets, similar to free-expanding Gaussian beam in laser beam physics. In laser beam physics, the 
spatial extent of the beam is well-defined by the location where the beam intensity drops to l/e 11 of the 
peak intensity, where n is usually 2. Connecting these locations as a function of axial distance for various 
values of n creates a form of ray tracing where all the rays lead back to the effective point of origin for the 
beam. While the method was interesting, difficulties related to analyzing the plume near the central axis 
of the thruster where the cathode plume dominates put some of the results into question. 

Having listed known near-field Faraday probe studies on the Hall thruster, this paper will now turn to 
how data analysis had been carried out in these studies. The focus in this section will be on how the point 
of origin and limits of integration are determined. 

The current density is a scalar value and the direction of the local ions cannot be determined from the 
current density. To calculate divergence angle from current density data, one must know where said 
current originates from. Of the mentioned prior studies that attempted to determine the divergence angle 
from near-field Faraday probe data, two studies (Refs. 8 and 13) assumed the point of origin is located 
along the centerline of the discharge channel at the channel exit plane. This assumption is invalid because 
interior probing has already revealed that the acceleration zone of the Hall thruster is located near the exit 
plane and has a non-negligible spatial extent from the perspective of the near-field region (Ref. 15). Or in 
other words, if one assumes the plasma originates from the intersection between the discharge channel 
centerline and the exit plane, then one cannot also assume that the source is a singular point. The only 
study that used plasma properties to determine the point of origin was only able to analyze the plasma 
occupying radial positions greater than the channel centerline (Ref. 1 1). As such, there is currently no 
proven method for finding the point of origin and is likely the reason why near-field Faraday probe data is 
generally not used to calculate plume divergence. Nevertheless, studying how the annular plasma jet from 
the Hall thruster transition from the near-field to the far-field requires an understanding of the direction 
and spread of said plasma jet. 

At this point in the paper, it is beneficial to discuss why near-field Faraday probe analysis generally 
requires some form of integration limits. Far-field Faraday probe study, the pre-cursor to near-field 
Faraday probe study, is usually carried out by sweeping the probe in a semi-circle around the thruster 
through the central axis of the thruster, also called a polar sweep. When swept this way, the surface 
normal vector for the collector of the probe is always pointed at the Hall thruster. Integration can then be 
carried out from -90° to +90° with the assumption that any ions collect beyond those angles are going the 
wrong way for thrust generation and likely a result of charge-exchange reaction. Further details on how to 
carry out a far-field Faraday probe study and deal with the effect of charge-exchange reaction can be 
found in Brown’s works (Refs. 16 and 17). A near-field Faraday probe cannot be swept in the same 
fashion as a far-field Faraday probe. Since most Hall thrusters have relatively flat faces, sweeping in the 
near-field is generally carried out either radially or axially in a rectangular coordinate system. 

Equation (1) shows the mathematical definition of the cosine of the divergence angle for a polarly-swept 
probe. 


? r /2 

271 R 2 j(0)cos0sin0d0 

<cos5> i= n (1) 

J 0 f7U/2 V J 

271 R 2 j(0)sin0d0 

Jo 

Where R is the distance from the thruster to the probe, j is the current density, and 0 is the polar angle 
defined as equal to 0 when pointed along the central axis of the thruster. Note the denominator of 
Equation (1) is equal to the integrated beam current. For a radially swept probe where the pointing 
direction of the probe is fixed, the equation for the cosine of the divergence angle becomes Equation (2). 
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Where j z has the subscript z to denote the fact that the probe is always pointed in the axial direction, r is 
the radial coordinate, z p is the axial position of the point of origin. Radial position of the point of origin is 
assumed to be 0. Polar sweeping of a near-field Faraday probe is not done because the probe cannot get as 
close to the thruster exit plane as in the case of rectangular sweeping and cannot be made to always face 
the local plasma (origin can no longer be treated as one point). Comparing Equation (2) to 
Equation (1) it becomes clear that while polarly-swept data is integrated from 0 to a finite number, 
radially-swept data is integrated from 0 to infinity. 

Aside from the mathematical issues associated with integration to infinity, two other experimental 
effects complicate the integration of radially- swept probe data. One, charge-exchanged (CEX) ions 
associated with a non-zero background pressure (separate from CEX ions associated with thruster 
neutrals) adds current to the measurement of a Faraday probe. At large values of r, the integration area is 
sufficiently large that even small CEX current can contribute to large total integrated current, even though 
this CEX current is not a part of the beam current. Reid’s work showed that integration of the current 
density from 0 to the limit of the data set often generate a total beam current that exceeded the discharge 
current, which is unphysical (Ref. 8). Two, noise and systematic offset in the data acquisition circuit can 
become greatly amplified when integrated to large values of radial position. To counter these undesirable 
effects, some kind of integration limits is generally necessary. 

Of the prior studies mentioned in this section, all but two studies used some scheme to determine the 
proper limits of integration when integrating for the total beam current. The two studies that are 
exceptions tried to subtract out the influence of the cathode plume, if applicable, and CEX ions, then 
integrated up to the limit of the measurement domain (Refs. 7 and 14). No detail is given on how the 
cathode plume is removed (Ref. 14). CEX ions were assumed to add a fixed amount of current density 
regardless of position and subtracted out. However, the production rate of CEX ions is known to depend 
on ion density, neutral density, ion energy, and distance from the thruster. The assumption of a constant 
CEX floor is questionable. In the remaining studies, the schemes for determining the limits of the 
integration can be divided into two categories, discharge channel integration and threshold-based limits 
(a.k.a. dynamic window integration). Discharge channel integration (Refs. 8 and 9) uses the radial 
positions of the inner and outer discharge channel walls as the lower and upper integration limits, 
respectively. Reid finds that discharge channel integration is only meaningful when the data is extremely 
close to the exit plane (Ref. 8), where plasma potential data indicates that some acceleration is still 
occurring (Ref. 15). Threshold-based limits (Refs. 8, 1 1, and 13) define the integration limits as the radial 
positions where the current density has dropped to some fraction of the peak current density. For traces 
with noticeable cathode plume, care must be taken during implementation to ensure the peak current 
density is that of the beam ions. For a free-expanding plasma jet, threshold-based limits will always 
bound the same amount of current density at every axial position and is the closest scheme to a physics- 
based scheme out of all those proposed for determining integration limits. Threshold-based limits also 
scales very well from one thruster to another as the absolute value of the threshold scales with peak 
current density, which is assumed to scale with discharge current. Nevertheless, the scheme of threshold- 
based limits has one major weakness; how one chooses the correct value of the fraction that determines 
the threshold is unclear. One possible approach is to calibrate the near-field Faraday probe results using 
far-field Faraday probe results. 

Proper choice of the integration limits is an ongoing area of research and is not resolved by the new 
analysis method. The authors believe more data is necessary to satisfactorily resolve this issue. 
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The New Method 


The following sections describe the various steps in the new near-field Faraday probe analysis 
method. To better illustrate the process, near-field Faraday probe data taken from the NASA-3 00M 
operating at 500 V discharge voltage and 20 kW discharge power is used as an example. Detailed 
description of the experimental setup and near-field Faraday probe testing of the NASA-300M can be 
found in a previous work (Ref. 11). All Hall thruster operating conditions are labeled as www V, yy kW, 
where www is the discharge voltage in volts and yy is the discharge power in kilowatt. All data presented 
are plotted on spatial coordinates that have been normalized by the mean channel diameter of the 
associated Hall thruster. R is the normalized radial coordinate and is 0 at the central axis of the thruster. R 
is positive traveling towards the 3 o’clock position of the thruster when the thruster is viewed face-on. Z 
is the normalized axial coordinate and is 0 at the thruster channel exit plane. Z is positive traveling 
downstream of the exit plane. Mean channel diameter is the average of the inner and outer discharge 
channel wall diameters. Thus, the centerlines of the discharge channels correspond to R = ±0.5. 

Note that for the purpose of this paper, only the near-field plasma upstream of plume merging is 
referred to as jets because they tend to be well-defined in spatial extent in the near-field close to the exit 
plane. 


Cathode Plume Removal 

Before applying the new near-field Faraday probe analysis method, the cathode plume needs to be 
treated. For data sets taken on Hall thrusters with centrally mounted cathodes, like the sets that will be 
shown later, the cathodes produced dense plasmas in the near-field that was recorded by the Faraday 
probe. Ions in the cathode plume are not beam ions and must be excluded from the divergence and beam 
current calculations. For the purpose of this study, the cathode plume is assumed to be an expanding jet 
that does not interact with itself or the beam plasma. However, the plume is allowed to be stretched 
slightly in the radial direction by any electric and magnetic fields present. The exact form that best fit the 
cathode plume is unknown so several different forms were tried. 

Since the cathode plume is known to interact with the beam plume at some downstream distance 
(Hall thruster operation require electrons from the cathode to mix with the beam plasma), the 
aforementioned assumption is only applicable for cathode plume very close to the exit plane. One must 
define a domain over which the cathode plume is curve-fitted. After performing analysis on the data from 
the 300M and 457Mv2 at many operating conditions, the integrated current from the region near the 
cathode was observed to rise as one traveled downstream. This trend is presumably due to increases in the 
mixing and interaction of the cathode plume with the beam plume. Another possibility is that the 
increased current is from a phenomenon known as the “central spike” observed in many other Hall 
thruster studies (including ones without a centrally-mounted cathode) (Refs. 7 and 8). At this point, the 
nature of this rise in integrated current is unclear so the fitting domain is limited to a relatively restrictive 
Z < 0.2. Within the limit of Z < 0.2, there are typically -10 radial sweeps. 

The radial boundaries for the curve-fit domain were chosen somewhat arbitrarily. As long as the 
domain captured most of the cathode plume but almost no beam plume, the curve-fitting results would be 
accurate. The radial boundaries are chosen to be a pair of lines radiating outward at ±25° from the orifice 
of the cathode plus a bit of extra space corresponding to twice the orifice diameter. Figure 1 shows the 
flooded current density contour for the 300M operating at 500 V, 20 kW. Overlaid on the figure are red 
dashed lines showing the boundaries of the cathode curve-fit region. 

Curve-fitting was performed on each of the Faraday probe sweeps that fell within the cathode curve- 
fit domain. Curve-fit forms that were tried include the Gaussian, the variable exponential, and the twin 
Gaussian. The authors quickly discovered that the Gaussian form left behind too much plume near the 
base of the cathode plume. The twin Gaussian form fitted well but used too many variables. The width of 
the twin Gaussian form was difficult to define since it involves two overlapping Gaussians. Same was 
also true of the variable exponential form. However, after performing curve-fits on the data from the 
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Figure 1. — Flooded current density contour with dashed 
boundary lines showing the cathode curve-fit domain. 
(300M operating at 500 V, 20 kW). 



R 


Figure 2. — Curve-fitting the Gaussian form to the 
cathode plume. (300M operating at 500 V, 20 kW). 
Black solid line is the raw trace, red dashed line is 
fitted curve, and blue solid line is the left-over when 
the red dashed line is subtracted from the black 
solid line. 


300M and 457Mv2 at many operating conditions, the authors noticed the power on the variable 
exponential form fell somewhere between 0.6 and 1.2 for all fitted traces with an average roughly around 
0.9. The analysis program is capable of handling a power of 0.9 but for the sake simplicity and to make 
analytic checks possible a power of 1 was used instead. Thus, the best fit form for the cathode plume in 
the region very close to the exit plane was found to be a simple exponential decay function. Equations (3) 
to (6) show the mathematical expressions for the various curve-fit forms used. 


Exponential decay: 

f(r) = A*exp[- r | / B] 

( 3 ) 

Gaussian: 

f (r) = A * exp[-r 2 / B] 

( 4 ) 


Variable Exponent: 

f (r) = A * exp[- r n / B] 

( 5 ) 

Twin Gaussian: 

f (r) = Aj * exp[-r 2 / ] + A 2 * exp[-r 2 / B 2 ] 

(6) 


Where A, B, n, Ai, Bi, A 2 , and B 2 are fitting constants. 

Figure 2 shows an example of curve-fitting the Gaussian form to the cathode plume. Note the bump 
near the base of the cathode plume in the left-over trace (blue solid line, grey when printed in black and 
white). 

Once curve-fit has been performed on the near- field plume close to the exit plane, the cathode plume 
was extrapolated downstream to the rest of the measurement domain. To perform this extrapolation, an 
overall curve-fit was found that adequately describe the individual curve-fit result. In particular, the width 
parameter, B, was plotted against the normalized axial coordinate, Z, and fitted to a second-order 
polynomial. For a free-expanding jet under no external force, the plot of B versus Z will form a straight 
line. By using a second-order polynomial, the jet was allowed to have some perturbation in the radial 
direction that could result from influence by external electric and magnetic forces. The integrated current 
of the cathode plume was also calculated and averaged over the curve-fit domain. Using the result of the 
polynomial fit of B versus Z and the average integrated cathode current, the fitting parameter A was re- 
calculated for the entire measurement domain. This step effectively re-normalized the equations fitted to 
the cathode plume so that the integrated cathode current was constant for all axial positions. The resulting 
overall fit follows the original assumption that the cathode plume is an expanding jet that does not interact 
with itself or the beam plasma. 
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Figure 3. — Examples of curve-fitting the exponential 
decay form to the cathode plume at various axial 
positions. (300M operating at 500 V, 20 kW). Black 
solid lines indicate the raw trace. Red dashed lines 
indicate the curve-fit results. CD stands for current 
density. 
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Figure 4. — The curve-fit parameters A, B, and the 
integrated cathode current as a function of axial 
position. Black dots indicate curve-fit results. Red 
dashed lines indicate the overall fit that will be used to 
project the cathode plume beyond the curve-fit domain. 



Figure 3 shows examples of curve-fitting the exponential decay form to the cathode plume at various 
axial positions. The black lines show the data and the red dashed lines show the curve-fit results. Note 
that sometimes the cathode current density was so high that the data acquisition system saturated (this 
occurred at -380 mA/cm 2 for the data set shown). When this occurred, the data analysis program 
excluded the saturated region while performing the curve-fit. Since the curve-fit result can be much 
higher than the raw data in the saturated region, subsequent subtraction of the curve-fit result from the 
raw data can cause artificial negative data to be introduced. The program resolves this issue by zeroing 
out any negative value within the saturated region after the subtraction is performed. 

Figure 4 shows the curve-fit parameters A, B, and the integrated cathode current as a function of axial 
position. The black dots are the curve-fit results within the cathode curve-fit domain (Z<0.2). The red 
dashed lines show the overall curve-fit that will be used to extrapolate the influence of the cathode 
throughout the measurement domain. The middle plot shows the second order polynomial fit of B versus 
Z. The lower plot shows the total cathode current at each Z position as well as the average over the entire 
cathode curve-fit domain. The upper plot shows the fit constant A from the curve-fit and from the 
normalization calculation. The upper plot is effectively a check on the overall curve-fit results. 
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Figure 5. — Curve-fitting the exponential decay form to the 
cathode plume. (300M operating at 500 V, 20 kW). 



Figure 6. — Illustration of a free-expanding jet emanating 
from a point in 2D axisymmetric space. 


Having found an overall curve-fit to the cathode plume, the cathode plume curve-fit was subtracted 
from the raw data, which sometimes generated odd looking spikes. This was often due to the slight 
misalignment between the raw data and the curve-fit results. To make sure these spurious spikes do not 
affect the analysis results, an additional section of code is added to exclude these spikes from the next 
step of the analysis. This section of code looks for the local minima between the spike and the rest of the 
plasma plume and limits the integration boundary appropriately. The effectiveness of this section of code 
can be seen in a later contour plot with the integration limits overlaid. Figure 5 shows an example of 
applying the overall cathode plume curve-fit at one axial distance. Note the line representing the left-over 
trace contains one spurious spike. 


The Iterative Pathfinding Method 

Since real experimental data is accompanied by noise, the problem of fitting the path of a plasma jet 
to said data is mathematically overdetermined. There is more than one solution to the problem and none 
of them fits perfectly. To find the best fit solution, a relatively simple and fast approach is to establish 
error metrics and minimize the error iteratively. The iterative pathfinding method starts with one 
simplifying assumption, which is easiest to understand from a 2D axisymmetric view (a cross-section 
view of the Hall thruster). In this view, the annular plasma jet show up as two jets, one on the +R side and 
one on the -R side, and are the mirror image of each other. The iterative pathfinding method assumes that 
each of the two plasma jets can be approximated as originating from a single point. The point becomes a 
ring in a 3D view. No assumption is made about the angular distribution of the current density (i.e., the 
current density does not have to have a Gaussian distribution). Figure 6 illustrates this assumption. The 
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plasma jets are assumed to expand freely beyond the channel exit plane at some angle, called the travel 
angle, and have some effective jet divergence angle. Within this paper, 0 p j stands for the jet travel angle 
on the +R side, 0 n j stands for the jet travel angle on the -R side, S p j stands for the divergence angle of the 
jet on the +R side, and 5 n j stands for the divergence angle of the jet on the -R side. Travel angles are 
defined as 0 when pointed axially and positive when leaning towards the +R direction. Thus, a jet on the 
positive side with a positive travel angle and a jet on the negative side with a negative travel angle are 
both pointing radially outward. 

For a free-expanding jet, both the travel angle and the divergence angle are constant as the jet 
expands but only if the calculation is carried out using the true point of origin. If one carries out the 
calculation using a different point as the origin, the travel and divergence angles will appear to vary with 
axial position. One can build an error metric out of the variations in the travel and divergence angles and 
try out different origin locations until the error is minimized. For example, if the guess origin location is 
higher up in the drawing than the red circle in Figure 6, the upstream part of the blue jet will appear to 
have a large positive travel angle while the travel angle of the downstream jet will appear to be small. One 
can then move the guess origin location downward until the travel angle calculated at different axial 
positions agree. On the other hand, if the guess origin location is somewhere along the solid black travel 
arrow but to the left or right of the red circle, the travel angles along the entire jet will appear to be the 
same, but the divergence angles will not. One can then move the guess origin location along the travel 
axis until the divergence angles calculated at different axial positions agree. 

In the actual implementation of the iterative pathfinding method there are two core steps. In step one, 
starting with a guess on the point of origin, the travel and divergence angles are calculated over various 
axial positions. Since the jet from the +R and -R sides are expected to eventually merge, the free- 
expanding jet assumption is only good in the near-field. In fact, a good way to define the near-field region 
is the axial extent over which there is no jet merging. Analysis of the 300M and 457Mv2 data at many 
operating conditions show that the jets can start merging as early as Z = 0.5 depending on the discharge 
voltage (lower discharge voltage generally trends with earlier merging). The authors also noticed the 
values of the angles and total integrated current tend not to be very stable in the very near- field, typically 
Z<0.2. This effect can be seen in later plots. The jet analysis zone (separate from the cathode curve-fit 
zone) is defined as 0.2<Z<0.5 based on the preceding observations. Equation (7) is the integral equation 
for calculating the travel angle. Integration limits were replaced with appropriate limits to be described in 
a later section. Mirroring these limits to the -R side gives the equation for 0 n j. 


271 j z (r)tan0rdr 
tan0 Di = ^ 

HI -oo 

2n Jo jzOOrdr 


tan 0 = - 


r - iv 


z-z r 


( 7 ) 


Where r p and z p are the radial and axial coordinates of the point of origin, respectively. Note that in 
theory, cosine can be used instead of tangent. However, the mathematical definition for cosine is such that 
positive small angle and negative small angle both yield positive cosine values. Thus, inverting cosine to 
obtain 0 in the present coordinate system can sometimes yield angle values that are in the wrong 
quadrant. Tangent does not have this issue. 

Equation (8) is the integral equation for calculating the divergence angle of each jet. Note that this is a 
measure of the spread in the plume around the travel angle, which is allowed to be non-zero in the current 
formulation. 
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Replacing the limits of integration with the appropriate values on the -R sign and replacing 0 pj with 0 nj 
will give the equation for the -R jet. 

In step two, the variations in the travel and divergence angles with axial distance is used to compute 
the error metrics. To simplify the calculations, the travel angle is assumed to be small. This assumption 
implies that the radial coordinate of the origin has a much greater influence on the travel angle than the 
axial coordinate, and that the axial coordinate has a much greater influence on the divergence angle than 
the radial coordinate. Pairing up the two controls with the two variables greatly simplifies the feedback 
equations. Equations (9) and (10) show the feedback equations for determining the next guess on the 
point of origin. 


r P ,k+i =r p,k“Pi si H 


r r r ^ 

r p,k-l ~ r p,k 

E e, P j,k-i -E e,pj,k 


^0,pj,k 


( 9 ) 


f 

z p,k+l = z p,k-p2sign 

v 


z p,k-l 

E 8,pj,k-1 


\ 


z p,k 

_E 6,pj,k y 


E S,pj,k 


( 10 ) 


Where the subscript k indicates the k-th iteration, Pi and p 2 are tuning parameters to control the speed of 
convergence, Ee p j k is the error metric on the travel angle, and Es^k is the error metric on the divergence 
angle. The function “sign” is equal to -1 or +1 depending on the sign of the content. If the sign functions 
are replaced with the contents of the sign functions, the above two equations would be the same as 
Newton’s method for root solving. However, the overdetermined nature of the problem seemed to give 
rise to convergence issue for Newton’s method. Instead, the sign function was used to maintain the ability 
of Newton’s method to determine the right direction to guess while the tuning parameters were used to 
control convergence speed. Note that the tuning parameters also control the accuracy of the guess. Due to 
the presence of noise in the data, the feedback algorithm tends to dance back and forth around the true 
point of origin. The amplitude of this numerical oscillation and the convergence speed of the algorithm 
both increase with p. Too high of a value of P can cause the algorithm to oscillate wildly and not 
converge. As such, Pi and p 2 were individually tuned to achieve adequately small numerical oscillation 
amplitudes while maintaining a decent convergence speed. Note that this simplified feedback algorithm 
will still converge even if the travel angle is not small but will take longer to converge than a more 
general feedback algorithm that links each error metric to both coordinates. On the hand, this simplified 
algorithm is numerically more stable than a more general algorithm. 

Equations (11) and (12) show the error metrics used for this study. 


E e.pj 


E 6,pj 




<ii> 


( 12 ) 


Where 0 pj - is the mean travel angle, S pj - is the mean divergence angle, and N is the number of axial 
positions over which the angles were calculated. Only data from the jet analysis zone (0.2<Z<0.5) were 
used in the calculations. If the 4 th power is replaced with the 2 nd power and the quartic root replaced with 
the square root, Equations (11) and (12) become root-mean-square of the difference between the angles 
and the mean angles. The 4 th power was used to increase the sensitivity of the error to difference between 
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the guess location and the true point of origin. The addition of noise to the data set creates many local 
minima in the vicinity of the global minimum that can destabilize the convergence process. Trials were 
performed with both the 2 nd and the 4 th power error metrics and the feedback algorithm was much more 
stable with the 4 th power error metrics. 

Figure 7 shows an example of the convergence plot for the iterative pathfinding method. The data set 
is from the 300M operating at 500 V, 20 kW. The black lines show the successive guesses in the 
coordinates of the +R jet while the red lines (grey if printed in black and white) show the guesses for the 
-R jet. For the -R jet, the negative of the radial coordinate is plotted so that the difference between the 
absolute value of the radial coordinates of the two jets is easier to see. Note that currently the program 
will run 500 trials regardless of whether convergence is detected. The main purpose of the program was to 
test out the concepts behind iterative pathfinding. The convergence algorithm can likely be further 
refined. Figure 8 shows the total beam current, travel angle, and divergence angle at the end of the 
analysis. The limit of integration used for this example was the threshold-based limits with a fraction of 
0.002 (i.e., 0.2 percent of peak current density). From Figure 8, one can see that the values of the various 
parameters can be quite different just upstream of Z = 0.2 than in the jet analysis zone. This trend 
appeared at many conditions for both the 300M and the 457Mv2. One possible explanation is that the jets 
are still accelerating upstream of Z = 0.2 since part of the acceleration zone is known to reside beyond the 
exit plane for these thrusters (Ref. 12). Then again, it is unclear how this acceleration could affect the 
Faraday probe results the way it does. 



Figure 7. — Example convergence plot for the iterative 
pathfinding method. (300M operating at 500 V, 
20 kW). The black lines represent the origin for the 
+R jet and the red lines represent the origin for the - 
R jet. Note that the negative of the R value for the - 
R jet is plotted instead of the actual value for 
convenience. 



Figure 8. — Example plot of total beam current, travel 
angle, and divergence angle versus axial position. 
(300M operating at 500 V, 20 kW). In all subplots, 
vertical lines indicate the boundaries of the jet 
analysis zone. For the angle plots, the black lines 
represent the origin for the +R jet and the red lines 
represent the origin for the -R jet. 
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To calculate the overall divergence angle as seen from far away, the momentum of each jet is 
assumed to be decompose-able into the momenta of two collimated beams. The beams each carry half of 
the momentum of the associated jet and radiates out at plus/minus the divergence angle with respect to the 
travel angle of the jet. This approach is represented mathematically in the Equation (13). 


e pj+ 8 pj + v 


' 8 pj + Pnj + 8 nj + Pnj - 8 nj 


(13) 


Limits of Integration 

As previously mentioned, choice of the limits of integration is a big issue when analyzing near-field 
Faraday probe data. Threshold-based limits were selected for this study because of its scalability. To re- 
iterate, limits for threshold-based limits are selected by find the radial locations where the current density 
has dropped below a certain fraction (called the threshold fraction) of the peak current density in the 
plasma beam. Several values of the threshold fraction were tried for this study. Part of the results will be 
shown in a sensitivity study within the Results section. The fraction values were picked to be as low as 
possible without generating unphysical results. However, while experimenting with the fraction values, 
the authors noticed there were values below which the feedback algorithm will not converge. One 
possible explanation is that the noise in the data, which can be a big fraction of the mean current density 
at large values of R, can overwhelm the feedback algorithm. If the integration limits are not sufficiently 
small, the noise recorded at large radial positions, amplified by the associated element area in the 
integration equation, can adversely affect the stability of the convergence process. Further investigation is 
required on this phenomenon. 

Curiously, the minimum fraction of peak current density needed to achieve convergence for the 300M 
data set is slightly lower than for the 457Mv2 data set; Perhaps suggesting that the 457Mv2 data set is 
relatively noisier. The minimum value of the fraction needed to achieve convergence for the 300M is 
-0.002 (i.e., 0.2 percent of the peak current density), while that for the 457Mv2 is -0.005. To make direct 
comparison between the two thrusters possible, results will be shown for limits based on both 0.2 and 
0.5 percent of peak current density. Figure 9 shows an example of the flooded current density contour 
after the cathode plume has been subtracted out. The red dashed lines indicate the limits of integration as 
functions of axial position. Data set comes from the 300M operating at 500 V, 20 kW. Integration limits 
are set to 0.2 percent of the peak current density. 



Z 

Figure 9. — Flooded current density contour with dashed 


lines showing the limits of integration. (300M operating 
at 500 V, 20 kW). Integration limits set to 0.2 percent of 
peak current density. 
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Results 

Experimental Data Sets 

To validate the iterative pathfinding method for analyzing near- field Faraday probe data, data was 
taken from a prior study on the NASA-3 00M and NASA-45 7Mv2 Hall thruster (Ref. 11). The thrusters 
are briefly described below. Detailed description of the thruster and experimental setup can be found in 
the prior study (Ref. 11). 

The NASA-3 00M and the NASA-45 7Mv2 are magnetic-layer Hall thrusters designed based on the 
principles outlined in Manzella’s dissertation work (Ref. 18). Both thrusters were originally developed for 
high-specific-impulse missions (Ref. 19). The two thrusters have lens-type magnetic field topology 
(Refs. 20 and 21). Each thruster ran with a centrally-mounted hollow cathode derived from the discharge 
cathode for NASA’s Evolutionary Xenon Thruster (NEXT) (Ref. 19). For the data used in this paper, the 
magnetic field settings of the thrusters were optimized to give roughly symmetric shapes while 
maximizing anode efficiency. The cathode flow fraction was maintained at 8 percent. 

The nominal operating condition of the 300M is 500 V, 20 kW. The data sets that will be shown for 
the 300M span seven operating conditions with discharge voltages between 200 to 500 V, discharge 
currents between 20 to 67 A, and discharge powers between 10 and 20 kW. The anode efficiency of this 
thruster varies from 57 to 73 percent over the tested operating conditions. The nominal operating 
condition of 457Mv2 is 500 V, 100 A. The data sets that will be shown for the 457Mv2 span 13 operating 
conditions with discharge voltages between 200 to 500 V, discharge currents between 40 to 100 A, and 
discharge powers between 10 to 50 kW. The anode efficiency of the 457Mv2 varies from 55 to 70 percent 
over the tested operating conditions. 

Additional data from a far-field Faraday probe study of the 300M is also included in order to compare 
the differences between the two types of Faraday probe studies (Ref. 22). The operating conditions of the 
300M studied during the far-field Faraday probe test matches those studied during the near-field Faraday 
probe test. The data from the two types of Faraday probe were taken using experimental setups that were 
identical in every way except for the diagnostics in use. 

Current Density Contour Before and After Cathode Plume Removal 

This section shows flooded current density contour plots before and after the cathode plume had been 
removed. Generally, the cathode plume removal results can be placed under one of three categories 
depending on how much plume remains along the thruster centerline in the near-field region after 
removal; they are “clean”, “plume”, and “anomalous”. The “clean” category refers to conditions for 
which no noticeable plasma remains (except the occasional spurious spike) in the region downstream of 
the cathode, prior to plume merging. The “plume” category refers to what appears to be a plasma plume 
somewhat downstream of the cathode that remains after the cathode plume removal process is applied. 
This is different from a spurious spike because the plume is wider and denser than what one would expect 
from curve-fit misalignment. The program is able to exclude this plume from the analysis domain the 
same way it excludes spurious spikes. The “anomalous” category is similar to the “plume” category 
except the plasma plume that remains downstream of the cathode after the cathode plume removal 
process is denser, and the program is not able to exclude the “anomalous” plasma from the beam analysis. 
One possible reason why the program could fail to exclude this anomalous plasma is that the data may be 
too noisy in that region. For the figures in this section, the “before” plots are overlaid with red dashed 
lines showing the cathode curve-fit zone while the “after” plots are overlaid with red dashed lines 
showing the integration limits corresponding to 0.2 percent of peak current density. 
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Figures 10 to 12 show the flooded current density contour plots before and after the cathode plume 
removal process. Figure 10 corresponds to the 300M operating at 500 V, 20 kW, and falls in the “clean” 
category. Figure 12 corresponds to the 457Mv2 operating at 250 V, 20 kW, and falls in the “plume” category. 
Figure 12 corresponds to the 300M operating at 300 V, 20 lcW, and falls in the “anomalous” category. 
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Figure 10. — Flooded contours of current density before 
operating at 500 V, 20 kW. 
l r *150 
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Figure 12. — Flooded contours of current density before and after cathode plume removal for the 300M 
operating at 300 V, 20 kW. 
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For the 300M, data sets for 6 of the 7 tested conditions fall in the “clean” category; 1 data set falls in 
the “anomalous” category. For the 457Mv2, data sets for 6 of the 13 tested conditions can be considered 
“clean”, 6 more falls under the “plume” category, and 1 falls in the “anomalous” category. Since the 
program is able to exclude the plume in the “plume” category data sets, there does not appear to be any 
issues with applying the beam analysis on 18 of the 20 tested thruster/conditions. For the two data sets 
that fall in the “anomalous” category, the uncertainty in the analysis results is likely to be higher, though 
the actual total current represented by the anomalous plume is no more than a few percent of the beam 
current due to the associated integration element area being small. It is currently unclear how best to 
estimate the uncertainties in the results of the iterative pathfinding method. Instead of calculating the 
increase in uncertainties, the data sets that fall under the “anomalous” category will simply be noted in the 
results. 

The presence of a plasma plume downstream of the cathode even after applying the cathode plume 
removal process implies that the probe collects more plasma there than can be accounted for by a free- 
expanding jet. However, what is really interesting about the observed results is that it varies greatly across 
operating conditions. For example, of the 7 cases from the 300M, only 1 had a noticeable central plume 
after cathode removal. Yet, for the 457Mv2, many cases displayed this unusual characteristic. Of the 
13 cases from the 457Mv2, all 5 cases that operated at 100 A discharge current displayed a sizable 
downstream cathode plume, while both cases that involved 80 A operation displayed a small downstream 
cathode plume. The remaining 6 cases involved operating at less than 80 A discharge current and 
displayed no noticeable downstream cathode plume after applying cathode plume removal. It would 
appear that the left-over plasma plume downstream of the cathode is linked to the discharge current 
(possible the current density) that the thruster is operating at. 

One possible explanation for the presence of left-over plasma plume downstream of the cathode is 
that the cathode removal process does not operate correctly for a large cathode current. For instance, if 
there is significant plasma generation downstream of the cathode, particularly beyond Z = 0.2, the extra 
plasma would not be picked up by the cathode plume removal algorithm. The amount of near-field 
plasma potential data in the cathode region on the 300M and 457Mv2 is somewhat limited (Refs. 10 and 
12). Based on this limited data, the area downstream of the cathode is observed to contain a long strip of 
negatively biased region. This region is radially thin (-0.1<R<0.1) and stretches out to about Z = 1. 
However, this region is only about 5 to 10 V lower than the surrounding plasma and the potential drop is 
arguably not high enough to ionize a meaningful amount of plasma. 

Another possible explanation is that the left-over plume is part of the “central spike” phenomenon 
that has been observed in many Hall thrusters (Refs. 7 and 8). While the exact nature of the central spike 
is unknown, the presence of plasma in the region just downstream of the thruster along the thruster 
centerline has been confirmed even when the cathode is not centrally mounted (Ref. 7). 

A third possible explanation is that the merging of beam plasma occurs much further upstream when 
the cathode flow is high. Recall that the cathode flow fraction was fixed at 8 percent. Since conditions 
with high discharge current corresponded to high anode mass flow rates, they also corresponded to high 
cathode mass flow rates. The increased mass flow may have enhanced the mobility of the plasma in the 
vicinity, allowing the plume to merge sooner after exiting the thruster than otherwise possible. 

Total Ion Current and Divergence Angle 

In order to obtain an accurate estimate of the total ion current, the effect of secondary electron 
emission (SEE) must be addressed. Secondary electrons born on a negatively biased probe will accelerate 
away from the probe. This effect adds extra current to the probe measurement that is indistinguishable 
from the collected ion current. Faraday probes are typically constructed from materials with low SEE 
yield like molybdenum and tungsten (Refs. 16, 25, and 26). The probe in cited data set uses molybdenum 
(Ref. 1 1). In the past, researchers assumed the effect of SEE can be ignored if such materials were used. 
While the singly-charged xenon-induced SEE yield for molybdenum and tungsten are indeed negligible 
(0.013 to 0.022), the doubly-charged xenon-induced SEE yield are roughly 10 times higher than the 
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TABLE L— SUMMARY OF SEE DATA FOR XENON ION 
BOMBARDMENT OF MOLYBDENUM (REFS. 23 AND 24) 


Bombarding particle 

SEE yield of molybdenum 

Xe + 

0.022 

Xe 2+ 

0.20 

Xe j+ 

0.70 


singly-charged yield, and the triply-charged SEE yield are roughly 35 times higher than the singly- 
charged yield (Refs. 23, 24, and 27) Table 1 summarizes the SEE yield values for multiply-charged xenon 
bombarding molybdenum. The singly-charged and doubly-charged xenon-induced yields are averages of 
the SEE yield data for ion energies in the range of 200 to 800 eV in Hagstrum’s 1956 work (Ref. 23). The 
triply-charged xenon-induced yield of molybdenum could not be found in the literature. The value in 
Table 1 is a projected value based on the similarity in yield between tungsten and molybdenum. The ratio 
of triply-charged induced yield to doubly-charged induced yield for tungsten is 3.5, so the yield for 
molybdenum is approximated to be 3.5 * 0.2, or 0.7 (Ref. 24). 

In a past work, Huang, et al., devised an approximate formula to estimate the effect of SEE based on 
current fractions of multiply-charged species (Ref. 22). This formula is reproduced in Equation (14). 


j _ — £fp — _ j pp * y 

1 + Yeff k 1 + Yk 


(14) 


Where J is the actual ion current density, J FP is the current density measured by a Faraday probe, y e ff is the 
effective SEE yield, y k is the SEE yield associated with bombardment by the k-th species, and f2 k is the 
current fraction of the k-th species. Using Equation (14), we can write a new equation to show the 
definition of the effective SEE yield for the Faraday probe, shown in Equation (15). 


Yeff = 



(15) 


Using the Equation (14) requires knowledge of the current fraction, which is typically obtained using a 
Wien filter (a.k.a. ExB probe). However, current fraction measurements were only obtained for the 300M, 
and not the 457Mv2. Since the main purpose of this study is to determine the effectiveness of the iterative 
pathfinding method, an approximation will be used for the effective SEE yield. From past measurements, 
the average effective SEE yield for the 300M over the range of operating conditions treated in this study 
was found to be 0.077±0.012 (Ref. 22). The 457Mv2 was designed based on the same principle as the 
300M, operates at very similar current density, and the data set being treated span very similar discharge 
voltages. Thus, by assuming the effective SEE yield is 0.077 for all data treated in this study, the random 
error in the total current rises by -1 percent while a 7.7 percent systematic error is removed. 

Table 2 lists the operating conditions and the angle results from applying the iterative pathfinding 
method with 0.2 and 0.5 percent threshold-based limits. An asterisk is added in front of the discharge 
voltage for each of the two operating conditions with “anomalous” plume that remains after cathode 
plume removal. A dash is shown if applying the new method with the corresponding threshold-based 
limits does not result in convergence within 500 iterations or converge to an improbable location. For 
example, when the feedback algorithm is barely stable, it will sometimes converge to a positive value of 
Z, which base on known Hall thruster physics is highly improbable. The data set for the 457Mv2 
operating at 500 V, 50 kW is a unique case in that the threshold fraction for integration limits must be 
1 .5+ percent in order to achieve convergence. It is unclear why this one case requires such a high 
threshold fraction in order for the feedback algorithm to work stably. 
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TABLE 2.— SUMMARY OF TRAVEL AND DIVERGENCE ANGLES FOR 
TWO DIFFERENT THRESHOLD-BASED LIMITS 


Thruster 

Disch. voltage, V 

Disch. current, A 

Disch. power, kW 

Overall 
divergence 
angle, deg. 

Travel angle, 
deg. 

(+R side) 

Divergence 
angle, deg. 
(+R side) 

Travel angle, 
deg. 

(-R side) 

Divergence 
angle, deg. 
(-R side) 

0.2% 

0.5% 

0.2% 

0.5% 

0.2% 

0.5% 

0.2% 

0.5% 

0.2% 

0.5% 


200 

50.3 

10.1 

21.6 

19.5 

8.0 

4.5 

22.0 

19.8 

-6.6 

-3.4 

21.2 

19.2 


300 

33.3 

10.0 

11.6 

8.6 

4.5 

2.7 

11.1 

8.5 

-3.5 

-1.2 

12.1 

8.7 


*301 

67.3 

20.3 

16.5 

15.2 

5.5 

3.8 

16.4 

15.2 

-4.3 

-2.8 

16.6 

15.2 

300M 

401 

24.9 

10.0 

12.9 

9.4 

4.7 

2.3 

12.3 

9.1 

-4.2 

-1.2 

13.5 

9.6 


399 

50.0 

20.0 

10.1 

8.5 

1.9 

1.0 

10.0 

8.5 

-0.9 

0.2 

10.2 

8.5 


501 

19.9 

10.0 

14.4 

10.4 

5.4 

2.4 

13.7 

10.1 

-5.2 

-1.3 

15.0 

10.6 


501 

40.6 

20.3 

12.0 

10.4 

3.4 

2.3 

12.2 

10.7 

-2.6 

-1.4 

11.8 

10.1 


202 

49.1 

9.9 

- 

13.6 

- 

5.3 

- 

13.4 

- 

-3.1 

- 

13.7 


*200 

99.8 

20.0 

- 

24.7 

- 

27.2 

- 

16.0 

- 

-22.3 

- 

14.8 


250 

80.2 

20.1 

- 

15.7 

- 

4.9 

- 

15.2 

- 

-4.3 

- 

16.1 


250 

99.8 

25.0 

- 

17.6 

- 

6.4 

- 

17.3 

- 

-6.4 

- 

17.9 


301 

66.4 

20.0 

- 

12.5 

- 

1.9 

- 

13.6 

- 

0.2 

- 

11.4 


300 

99.8 

29.9 

- 

16.1 

- 

5.5 

- 

17.4 

- 

-1.9 

- 

14.8 

457Mv2 

400 

49.7 

19.9 

- 

11.4 

- 

5.7 

- 

12.0 

- 

-2.3 

- 

10.8 


400 

74.8 

29.9 

- 

12.5 

- 

-0.4 

- 

13.3 

- 

1.5 

- 

11.7 


401 

98.2 

39.4 

- 

10.9 

- 

3.3 

- 

10.1 

- 

-3.1 

- 

11.8 


501 

40.5 

20.3 

- 

12.7 

- 

6.4 

- 

13.4 

- 

-2.4 

- 

12.0 


500 

60.1 

30.1 

- 

10.6 

- 

4.1 

- 

10.6 

- 

-1.9 

- 

10.6 


501 

79.7 

39.9 

- 

11.7 

- 

3.4 

- 

11.6 

- 

-1.6 

- 

11.9 


501 

101.1 

50.7 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 


Studying Table 2 yields some insights into the new method. Using smaller threshold fraction for the 
integration limits correspond to larger angle values. This is expected because smaller threshold fraction 
means current density at larger radial position will be used in the angle calculations. The travel angle is on 
the average within a few degrees (except for one case) and is positive on the +R side, negative on the -R 
side. This trend implies that there is a small outward bias in the travel angles of the plasma jets. One 
possible explanation for the small outward bias is the swirl in the ions applied by the magnetic field. 
Alternately, the outward bias in the travel angle could be a systematic error in the new analysis method. 
Either way, the travel angle is small enough that associated cosine thrust loses are negligible. The case of 
457Mv2 operating at 200 V, 20 kW shows a noticeably higher travel angle than the other cases. Trying 
several different threshold fraction values showed that the program is indeed converging properly. It is 
unclear at this point whether the large travel angles are real. 

Table 3 lists the total beam current, current utilization efficiency, overall divergence angle, and 
divergence efficiency by operating conditions. The current utilization efficiency and the divergence 
efficiency are defined in Equations (16) and (17), respectively. The current utilization efficiency describes 
loss in efficiency associated with part of the discharge current being carried by electrons instead of ions. 
The divergence efficiency describes loss in efficiency associated with plume divergence. Both are 
typically calculated from Faraday probe data. 
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Where lb is the total beam current and Id is the discharge current. From Table 3 we see that using smaller 
threshold fraction for the integration limits correspond to larger total beam current, likely due to an 
expansion of the integration domain. For the lack of a better method to determine the right threshold 
fraction value, this study will rely on far-field Faraday probe results to calibrate the near-field Faraday 
probe results. For state-of-the-art Hall thrusters, this percentage is typically 70 to 80 percent, (Refs. 8, 28 
to 30) so a fraction value of 0.2 percent is a closer fit than 0.5 percent. However, cases for the 457Mv2 
generally did not converge when a fraction value of 0.2 percent was used. Further investigation with a 
refined convergence algorithm and more comparisons between near-field and far-field Faraday probe data 
is needed to determine whether 0.2 percent is the correct threshold value to use. 


TABLE 3.— SUMMARY OF BEAM CURRENT AND DIVERGENCE ANGLE FOR 
TWO DIFFERENT THRESHOLD-BASED LIMITS 


Thruster 

Disch. voltage, V 

Disch. current, A 

Disch. power, kW 

Total beam 
current, A 
(Corrected for 
SEE) 

*lb 

Overall 
divergence 
angle, deg. 

*ld 

0.2% 

0.5% 

0.2% 

0.5% 

0.2% 

0.5% 

0.2% 

0.5% 


200 

50.3 

10.1 

40.5 

31.7 

0.750 

0.630 

21.6 

19.5 

0.864 

0.889 


300 

33.3 

10.0 

26.4 

22.1 

0.737 

0.664 

11.6 

8.6 

0.960 

0.978 


*301 

67.3 

20.3 

53.9 

43.6 

0.745 

0.648 

16.5 

15.2 

0.920 

0.931 

300M 

401 

24.9 

10.0 

19.5 

16.2 

0.728 

0.652 

12.9 

9.4 

0.950 

0.974 


399 

50.0 

20.0 

40.3 

33.9 

0.750 

0.678 

10.1 

8.5 

0.969 

0.978 


501 

19.9 

10.0 

16.3 

13.5 

0.763 

0.677 

14.4 

10.4 

0.938 

0.968 


501 

40.6 

20.3 

32.8 

28.1 

0.752 

0.692 

12.0 

10.4 

0.957 

0.968 


202 

49.1 

9.9 

- 

30.9 

- 

0.630 

- 

13.6 

- 

0.945 


*200 

99.8 

20.0 

- 

69.5 

- 

0.696 

- 

24.7 

- 

0.825 


250 

80.2 

20.1 

- 

52.1 

- 

0.649 

- 

15.7 

- 

0.927 


250 

99.8 

25.0 

- 

65.8 

- 

0.659 

- 

17.6 

- 

0.909 


301 

66.4 

20.0 

- 

41.7 

- 

0.628 

- 

12.5 

- 

0.953 


300 

99.8 

29.9 

- 

67.0 

- 

0.671 

- 

16.1 

- 

0.923 

457Mv2 

400 

49.7 

19.9 

- 

33.0 

- 

0.664 

- 

11.4 

- 

0.961 


400 

74.8 

29.9 

- 

50.4 

- 

0.674 

- 

12.5 

- 

0.953 


401 

98.2 

39.4 

- 

59.6 

- 

0.607 

- 

10.9 

- 

0.964 


501 

40.5 

20.3 

- 

28.8 

- 

0.710 

- 

12.7 

- 

0.952 


500 

60.1 

30.1 

- 

39.9 

- 

0.665 

- 

10.6 

- 

0.966 


501 

79.7 

39.9 

- 

50.3 

- 

0.631 

- 

11.7 

- 

0.959 


501 

101.1 

50.7 

- 

74.5 

- 

0.737 

- 

- 

- 

- 
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Figure 13. — Sensitivity study of the divergence angle 
versus the threshold fraction for several operating 
conditions of the 300M. 


Self-Consistency of the Iterative Pathfinding Method 

To determine the validity of the iterative pathfinding method, this study will look at the self- 
consistency of the method and comparison to far-field Faraday probe results, which is in the next sub- 
section. The sensitivity study was performed by varying the limits of integration over a range of values 
and seeing how the analysis method responded. As previously mentioned, the current convergence 
algorithm became unstable at very low values of the threshold fraction. Non-convergent results will not 
be shown. Since the iterative pathfinding method was devised to tackle the issue of calculating point of 
origin and divergence angle, the divergence angle will be the primary metric of self-consistency. 

Figure 13 shows the divergence angle calculated by the iterative pathfinding method as a function of 
threshold fraction for four cases of the 300M. All seven cases displayed similar behavior but some were 
not plotted to avoid clutter. In general, the divergence angle was found to be within ±2° of the mean over 
a fairly large range of the threshold fraction. The ability of the iterative pathfinding method to determine 
angles appeared to be self-consistent. 

Comparison With Farfield Faraday Probe Data 

Comparison between near-field and far-field Faraday probe data is carried out using data sets from 
the NASA-300M. The far-field Faraday probe data were obtained using polar sweeps at a distance of 
5 mean channel diameters away from the thruster. CEX ions in the wings of the sweeps were removed by 
performing exponential curve-fits from 10° out to 25° and then extrapolating to 90° (Ref. 22). SEE effects 
had been corrected from the results. 

Table 4 compares the results of near-field Faraday probe analysis to far-field Faraday probe analysis. 
The threshold fraction for this set of near-field Faraday probe results is 0.2 percent. While performing the 
sensitivity analysis, the authors note that a threshold fraction of 0.2 percent also yield a very good 
agreement on the total beam current, in addition to being the lowest fraction for which applying the 
iterative pathfinding method on the 300M data would converge. This trend is likely a coincidence as it 
does not hold true for the 457Mv2 data. From Table 4 one can see the good agreement in total ion current 
between using the near-field Faraday probe data and the far-field Faraday probe data, which leads to very 
similar current utilization efficiencies. The calculated divergence angles are generally lower for the near- 
field Faraday probe data, which leads to higher divergence efficiencies. As result, the products of the two 
efficiency factors are higher for the near-field results but only by ~ 5 percent. The results in Table 4 
suggests that applying iterative pathfinding method to near-field Faraday probe data is generating the 
same trends and very similar values as the far-field Faraday probe analysis method. 
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TABLE 4.— COMPARISON BETWEEN NEAR-FIELD AND FAR-FIELD FARADAY 
PROBE ANALYSIS FOR THE NASA-300M 


Disch. voltage, V 

Disch. current, A 

Disch. power, kW 

Total beam 
current, A 
(Corrected 
for SEE) 

rib 

Divergence 
angle, deg. 

rid 

T|b*r|d 

NF 

FF 

NF 

FF 

NF 

FF 

NF 

FF 

NF 

FF 

200 

50.3 

10.1 

40.5 

40.1 

0.806 

0.806 

21.6 

26.3 

0.864 

0.804 

0.697 

0.648 

300 

33.3 

10.0 

26.4 

25.7 

0.793 

0.766 

11.6 

15.1 

0.960 

0.932 

0.761 

0.714 

301 

67.3 

20.3 

53.9 

52.6 

0.801 

0.776 

16.5 

20.4 

0.920 

0.878 

0.736 

0.681 

401 

24.9 

10.0 

19.5 

19.7 

0.783 

0.777 

12.9 

14.9 

0.950 

0.934 

0.744 

0.726 

399 

50.0 

20.0 

40.3 

38.2 

0.807 

0.768 

10.1 

14.4 

0.969 

0.938 

0.782 

0.720 

501 

19.9 

10.0 

16.3 

16.1 

0.820 

0.795 

14.4 

15.3 

0.938 

0.930 

0.769 

0.740 

501 

40.6 

20.3 

32.8 

31.4 

0.808 

0.780 

12.0 

14.7 

0.957 

0.936 

0.773 

0.729 


One way to interpret the similarity between the near-field and far-field Faraday probe results is that 
the far-field plume is simply a super-positioning of the near-field plasma jets. This interpretation is in 
agreement with the mean-free-path calculated from the plasma density in the transition region (0.5<Z<5). 
On the other hand, one can also interpret the small difference between the near-field and far-field Faraday 
probe results to be a sign of plume-merging action. If one take the jets from the two sides of the thruster 
and intersect them, one would expect some level of scattering, which increases plume divergence. On the 
average, the plume divergence angle found in the far-field was 3° higher than in the near-field, which may 
reflect the scattering effect. 

The results in Table 4 also bring up a previously mentioned issue for near-field Faraday probe 
analysis. While the total beam current calculated from near-field and far-field data agrees well, the values 
can be tweaked by adjusting the threshold fraction. For example, decreasing the threshold fraction 
(assuming convergence is still achieved) would extend the integration domain, increase the total beam 
current, and increase the divergence angle. The problem is that there is no physical justification for 
selecting any given value of the threshold fraction other than to match far-field Faraday probe data. Thus, 
the value of a near-field Faraday probe for engineering analysis of a Hall thruster is unclear when a far- 
field Faraday probe gives less ambiguous results and is generally easier to implement. However, near- 
field Faraday probe was originally designed for scientific studies of the evolution of the near-field plume 
and appears well suited for that purpose. The near-field Faraday probe is able to map out Hall thruster 
plasma structures that are otherwise undetectable. The observation earlier in the paper on “anomalous” 
plasma plume downstream of the cathode is one example of how the near-field Faraday probe can help 
Hall thruster researchers discover new phenomena in the near-field region. 


Summary and Future Works 

A new analysis method, called iterative pathfinding, had been devised for analyzing near-field 
Faraday probe data. Unlike most previous methods, the iterative pathfinding method relies on the 
properties of the plasma to determine the location of the point of origin, from which travel and divergence 
angles of the plasma jet can be calculated. The iterative pathfinding method was applied, after removal of 
the cathode plume, with threshold-based integration limits. Sensitivity analysis of the divergence angle 
results over a range of threshold fraction showed that method is self-consistent. Comparison of the near- 
field Faraday probe results to far-field Faraday probe results for the 300M showed good agreement. The 
total beam current calculated from the two sets of data agreed when a threshold fraction of 0.002 is 
selected. At the same threshold fraction, the divergence angle found in the near-field was on the average 
3° lower than in the far-field. The difference in divergence angle could be a result of ion scattering in the 
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transition region between near- field and far- field. However, the authors noted that the choice of 
integration limit can affect the near-field Faraday probe results, particularly the total beam current. 

A method for removing the plume of a centrally-mounted cathode from the near- field Faraday probe 
data was also tested. A dense plasma plume was observed downstream of the cathode along the thruster 
central axis even after applying cathode plume removal for operating conditions with high discharge 
current. Several possible explanations were proposed though none could be confirmed or denied in this 
study. 

The next step to improve the iterative pathfinding method is to refine the convergence algorithm. The 
goal is to improve the method to the point where it can be applied to various data sets, even ones with a 
lot noise. Another important next-step is to; if possible, devise a mean of overcoming the ad hoc way in 
which integration limits are set. 
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